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INITIATION OF COOLING DUE TO BUBBLE 


GROWTH ON A HEATING SURFACE* 
by Robert C. Hendricks and Robert R. Sharp 
Lewis Research Center 


SUMMARY 

The surface temperature profile of a nickel -chrome -alloy heater strip 
during nucleate boiling was matched to high-speed photographic data. The ini- 
tiation of surface cooling due to bubble growth was investigated and was dis- 
cussed in relation to existing models of nucleate boiling. Experimental data 
indicate that, at moderate subcooling, evaporation from a surface liquid film 
accounts for most of the ebullient heat transfer. No increase in heat transfer 
due to microconvection or surface quenching by cold bulk: liquid was observed. 


INTRODUCTION 

Much of the research concerning nucleate boiling during the' past several 
years has been conducted in order to determine the mechanisms, evidently due to 
bubble formation, that are responsible for producing the high heat fluxes 
associated with boiling. A quantitative description of the heat-transfer pro- 
cess has by no means been accurately formulated at present. Of the several 
models of boiling that have been proposed, the most widely accepted models 
hypothesize that the bubble acts as an agitator of the liquid. 

Rohsenow and Clark (ref. l), Gunther and Kreith (ref. 2), and others have 
suggested that the principal path of heat transfer is from the surface through 
the turbulent liquid between bubbles. Since conduction through a laminar 
thermal layer is insufficient to account for observed heat fluxes, these 
workers have proposed that growing and collapsing bubbles create random micro - 
convection currents of high velocity along the heated wall. In effect, agita- 
tion from the bubble is believed to destroy the thermal layer beneath the 
bubble and over some unknown distance around its base; this area of effective 
agitation is commonly referred to as the area of influence of a bubble. In 
the work of reference 3, in fact, it was observed by means of schlieren motion 
pictures that the thermal-layer thickness is greatly diminished during bubble 
growth within a distance of one bubble diameter around the nucleation site. 


*An abbreviated version of this report was presented at the American 
Institute of Chemical Engineers meeting, Los Angeles, California, February 4-7, 
1962 . 



These experiments give a clear definition of the region of effective forced 
convection due to a bubble. 


Forster and Greif (ref. 4) have suggested that a more efficient heat- 
transfer mechanism occurs -when a growing bubble displaces the superheated 
liquid layer entirely from the surface (fig. 1(a)). Cooling then occurs dur- 
ing collapse when cold liquid suddenly quenches the hot surface. Thus the 
heat-transfer mechanism would be a continuous cycle of vapor-liquid displace- 
ment with the surface heat flux increasing locally during the quenching 
process. 



Hot liquid being pushed Surface cooling 

into bulk with bubble collapse 


(a) Surface cooling by agitation. Model of reference 4. 



(b) Surface cooling by continuous liquid-film 
evaporation with possible surface drying area. 
Dependent on surface conditions. Model of 
references 5, 6, and 10. 


Figure 1. - Surface cooling effects in relation to bubble models. 


Evaporation into the bubble is generally thought to occur from the liquid 
superheat layer pushed away from the surface by a growing bubble. As evapora- 
tion takes place, the layer cools and the rate of this process controls bubble 
dynamics. More recently, Snyder (ref. 5) and Edwards (ref. 6) have suggested 
that evaporation occurs chiefly from the base of the bubble, either from the 
perimeter of the base or from a thin liquid film left on the surface during 
bubble formation. Rapid evaporation from the base would have the effect of 
drawing large amounts of heat directly out of the boiling surface during bubble 
growth (fig. 1(b)). 


The contribution to the total heat transfer due to evaporation into the 
bubble has, in general, been considered small if not negligible. Rohsenow and 
Clark (ref. l) have made measurements of bubble volumes from motion pictures 
and have compared the latent heat associated with the observed vapor flux with 
the total heat flux. Their results indicate that evaporation accounts for no 
more than 1 to 2 percent of the total heat transfer. These workers, however, 
ran boiling tests in the subcooled regime and did not correct for condensation 
in the growing bubble. Bankoff (ref. 7) measured condensation rates of steam 
bubbles and concluded that the vapor evolved is probably 10 times greater than 
the amount observed by Rchsenow and Clark. This still indicates only a minor 
role for an evaporative mechanism in the total heat transport. 

Most important to the support of a bubble agitation theory, however, 
would be a direct indication that bubble dynamics is sufficient to increase 
heat fluxes to levels observed in nucleate boiling. It was reasoned in the 
study of reference 8 that streams of air bubbles that were electrolytically 
generated on a heated surface should simulate the agitation over a boiling 
surface while the evaporative heat transport is virtually eliminated. The air 
bubbles were indeed found to have a strong effect on the heat-transfer coef- 
ficient, but no heat fluxes comparable to vigorous nucleate boiling were ever 
attained. 
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Another , but quite indirect , source of support for an agitation theory is 
the application of several semi theoretical correlations (ref. 9) that are based 
on a forced convection (modified Nusselt number) relation. This type of corre- 
lation has given fairly good results . 

Moore and Mesler (ref. 10) sought more direct evidence concerning the 
mechanism of boiling heat transfer by observing the temperature -time profile of 
a boiling surface with a high-response thermocouple read by an oscilloscope. 
They noted regularly an exceedingly rapid drop in temperature followed by a 
smaller rapid rise and a gradual reattainment of the initial temperature. 

Moore and Mesler did not know whether these drops occurred with the same fre- 
quency as bubble formation, nor did they have any indication of where, with 
respect to bubble growth, the temperature drops began. The rapidity of cooling 
immediately suggested an evaporation phenomenon, although similar temperature 
behavior would be expected on the flooding of a hot surface with cold liquid. 
They calculated the maximum heat flux obtainable from surface flooding and 
found this to be insufficient to account for the magnitude of the temperature 
drop. They also calculated that between 70 and 90 percent of the total heat 
transfer could be attributed to the temperature drop. This indicates that, if 
the drop is indeed caused by evaporation, latent heat transport is many times 
as important in the boiling phenomenon as had been thought previously. Cer- 
tainly such a drop is an outstanding characteristic of the microlayer model of 
evaporation described by Snyder (ref. 5) and Edwards (ref. 6). 

Both the quantitative and the qualitative conclusions of Moore and Mesler 
are open to considerable doubt, however, since the construction of their 
thermocouple caused a heating discontinuity Still, it is evident that a 
rapid drop in surface temperature does occur during boiling and that this drop 
is probably associated with an appreciable portion of the total heat flux. 

In this report high-speed photographic data have been matched with the 
surface temperature-time profiles underneath large discrete bubbles in order 
to determine when, with respect to bubble life, local fluctuations of surface 
heat flux occur. In particular, it was desired to determine exactly where in 
the ebullition cycle the large temperature drop noted by Moore and Mesler 
begins. If the drop occurs at the incipience of bubble growth rather than 
during departure, it will be caused by evaporation at the bubble base. The 
agitation and the liquid microlayer hypotheses are discussed in relation to 
the observed temperature fluctuations during the ebullition cycle. Absolute 
magnitudes of the temperature drop were also measured, and the resultant local 


^The thermocouple used in the investigation of reference' 10 consists of an 
electrically insulated Alumel wire core passing through a Chromel cylinder, 
which in turn is attached to an electrically heated surface. While the Chromel 
cylinder dissipates heat at a rate comparable to the homogeneous heater, the 
Alumel wire core does not dissipate' electrical energy. Hence, the core wire is 
heated only by conduction from the surrounding material. Heat removal from 
this core wire during a temperature fluctuation would occur at a greater rate 
than from the Chromel cylinder, and a temperature gradient across the thermo- 
couple junction would result. Consequently, the thermocouple would not read 
correctly, although the thermocouple would respond well to temperature pulses. 
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heat fluxes were calculated. From these data several conclusions have been 
reached concerning the relative quantitative importance of the agitation and 
the evaporative heat -transfer mechanisms. 

The results of this experiment are applicable to the boiling cases of 
discrete bubbles and bubble columns when subcooling does not exceed about 
30° F. Conclusions regarding microconvection and quenching might be inaccurate 
at higher subcoolings-. The poorly understood regime of vapor columns also 
requires further investigation. 


SYMBOLS 


cs camera speed 

L heater length 

l heater thickness 

q heat flux 

q ge n generated heat flux 

T temperature 

T max maximum temperature 

T gat saturation temperature 

t time 

^rm 1 reference time for mechanical system with bubble occurring before 
reference mark 

t^ g reference time for mechanical system with bubble occurring after 
J reference mark 

t r o l reference time for optical system with bubble occurring before 
reference mark 

t ro 2 reference time for optical system with bubble occurring after reference 

9 mark 

t<p phase lag 

to system time response 

V voltage 

x number of frames from reference timing mark to initiation of bubble 
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y number of frames until bubble crosses thermocouple location 
z heater coordinate 
a thermal diffusivity 
6 temperature 

cp transformer phase change, sec 


APPARATUS MB PROCEDURE 



The data-recording system is shown schematically in figure 2. 
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Figure 2. - Data-recording system. 


In order to 
fluid was 
free of foreign materials, the 
vessel was purged of air and was 
loaded with distilled, degassed 
water. The desired degree of 
sub cooling was obtained by draw- 
ing a vacuum in the vessel. Two 
types of electrically heated 
Nichrome elements were used. The 
first was 0.0056 inch thick and 
l/l6 inch wide, and the second 
was 0.001 inch thick and l/8 inch 
wide. The heated length of each 
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element was 1— inches . The 
16 



Figure 3. - Thermocouple installations. 


5 


L 




Nichrome ribbons were mounted on a block that insulated the lower surface. 

The first element (0.0056 by — by 1— in.) was relatively heavy in con- 

16 16 

struct ion and was used to observe qualitatively the temperature fluctuations 
during the ebullition cycle. Three 36-gage copper-constantan thermocouples 
were silver soldered to the heater subsurface, and the excess silver solder 
was removed from the area indicated in figure 3(a). This method of attachment 
provides a surface with no discontinuity, but involves a heater nonhomogeneity . 
The absolute values of the temperature readings at the subsurface are not con- 
sidered reliable, but the subsurface readings corrected for time response are 
reliable in showing temperature fluctuations. 

i ’-j 

The second element (0.001 by — by 1-^- in.) was capable of measuring tem- 

8 16 

perature profiles with far greater accuracy than the first. Since the thermo- 
couple acted as a heat sink and disturbed the heat-flux pattern out of the 
strip, the 36-gage ( 0.005-in. -diam.) copper-constantan leads were replaced with 
0.0005-inch-diameter Chromel-Alumel leads. Also, since any excess silver 
solder formed a path of low resistance to the heater voltage and created an 
area of lower heat generation, the Chromel-Alumel thermocouple was spotwelded 
at a point to the heater strip. The junction was made by crossing a Chromel 
and an Alumel wire and by spotwelding them to the strip at their point of 
contact as is shown in figure 3(b). A drop of epoxy cement was placed on the 
weld to give the two thermocouples strength and also to keep the lead wires 
away from the heater strip. The thermocouple junction shown in figure 3(b) was 
localized to an area of approximately 0.001 to 0.002 inch in diameter. A 
photomicrograph of the thermocouple junction is shown in figure 4. Direct - 
current pickup from the heater voltage was rather large, even across this tiny 
junction, but was determined by suitable calibrations. 



0.001 in. 


(a) Light field illumination. 



C -63731 


(b) Dark field illumination. 


Figure 4. - Welded thermocouple junction. 


6 


During each run, heater voltage and current, pool temperature, and pres- 
sure were measured. Subsurface temperatures were recorded on an oscillograph, 
and motion pictures were taken simultaneously with a high-speed camera. 
Subatmospheric pressures yielded larger bubbles, which facilitated the study of 
single bubbles . 


Determination of Time Base 

Correlation of the time base is a problem that arises when data are 
recorded simultaneously by camera and oscillograph or oscilloscope. Time 
shifts may occur because of phase changes, recording-circuit lags, lamp re- 
sponse to applied voltage, heater thickness, or thermocouple current output. 


Time Lags in Sensing Temperatures 

When a temperature fluctuation occurs on the upper surface of the heating 
strip, a finite time passes before this fluctuation influences the subsurface 
temperature. An approximation to the subsurface temperature pulse may be 
obtained by reducing the general transient heat -conduction equation to 


where 


Se c> 2 e 


G = T z = l 



z = 0 


(i) 


e = o z = o t = o 

This equation does not include a heat -generation term but is considered ade- 
quate since the local heat flux during the rapid temperature drop is more than 
one order of magnitude greater than the average electrical heat flux (fig. 7). 

Solutions of equation (l) along with a parametric plot of T/t against z/l 

f p max * 

for successive values of cxt/l c are shown in figure 11 of chapter 3 in ref- 
erence 11. For a sudden change in temperature on the boiling surface, the 
ratio of subsurface (thermocouple) temperature to the maximum change in surface 
temperature may be determined at any given time. 

To sense a change in temperature, rather than its absolute value, the 
plots of reference 11 can be consulted to determine the time interval between 
the initiation of a pulse and the time when a measurable temperature change 
occurs on the subsurface. This time lag is approximately 1.1 milliseconds for 
the 0 .0056 -inch -thick heater, an appreciable correction, but only 0.036 milli- 
second for the 0.001 -inch -thick heater. 
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Since the thin heater was used to measure absolute values of the tempera- 
ture profile, it is important to show that, during any actual temperature fluc- 
tuation on the boiling surface, the subsurface temperature will approach the 
supersurface temperature in a very short time compared with the total duration 
of the temperature change. This temperature lag across the heater will have 
the greatest magnitude when an instantaneous temperature pulse is applied to 
the upper surface. The solution of this problem also can be found in the 
curves of figure 11 of chapter 3 in reference 11. For a pulse of any magni- 
tude, 90 percent of the total temperature change is observed on the subsurface 
within half a millisecond. By comparison with the temperature profile (fig. 7), 
which has a duration of 15 milliseconds, it is apparent that the supersurface 
drop is far from instantaneous and that this drop must be quite close to the 
one actually sensed by the thermocouple. 

The heat capacity of the thermocouple itself is not considered to be an 
appreciable source of error since the dimensions of the junction are small 
compared with the heater thickness. With certain corrections for the effects 
of electrical pickup, the temperature profiles sensed by the thermocouple are 
considered to yield accurate absolute values. 


Oscillograph Recording System 

During the course of the experiment three types of thermocouple recording 
systems were used. In the first system a 400-microfarad damping capacitor, 
which increased the system delay time but reduced pickup, was wired in parallel 
with a 40-cycle galvanometer. In the second system pickup was reduced and the 
large capacitor was eliminated; consequently, the time -delay error was lowered. 
In the third system, which was used for measuring absolute values of the tem- 
perature profiles, a direct -current differential amplifier impressed the ther- 
mocouple signal across a 400 -cycle galvanometer. Control experiments^ showed 
that this galvanometer was capable of much faster response than was required 
for any of the temperature fluctuations studied. The oscillograph, which was 
run at paper speeds of 50 inches per second, could be read easily to 0.01 mil- 
livolt and with relative accuracy to 0.005 millivolt. All known phase lags 
were eliminated where possible. 


Optical System 

The time lag of the film to the frame of light is very small; however, the 
reference timing light response depends on the applied voltage. When a 
60-cycle voltage is used, the lamp may lag the applied voltage by 60°. After 
several runs were taken with this system, the 60-cycle timer was replaced by a 
square-wave timer because of reading and phase uncertainties in the 60-cycle 
system. 


p 

In one experiment, a few drops of cold water were dropped on the heater 
surface at the thermocouple junction. The action was recorded by the high- 
speed camera, and the response was recorded by the oscillograph. The result 
was nearly a square wave response in temperature when the drop contacted the 
heater surface. In a second experiment, the response of a galvanometer to a 
square-wave audio-oscillator input was observed. 
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(b) Reference time between optical - and oscillograph-recording 
systems with bubble occuring after reference timer mark 
(fig. 5(a)). 


(c) Reference time between optical- and oscillograph-recording 
systems with bubble occuring before reference timer mark. 


Figure 5. - Time base. 


9 


L 



TABLE I 


OPTICAL AND MECHANICAL TIME CORRELATION MEASUREMENT 


Number 

Camera 

Number 

System 

Phase 

Other 

Optical 

Meehan- 

Bub- 

of 

speed. 

of 

time 

lag. 

phase 

system 

ical 

bles 

frames 

cs. 

frames 

response. 

• 

fc 

cpj- 

ec 

changes , 

reference 

system 


from 

frames/sec 

until 


ms 


9, 

time. 

ref- 


ref- 


bubble 

msec 




msec 

^ro ,2* 

erence 


erence 


crosses 








msec 

time. 


timing 


thermo- 









’ t rm,2 ' 
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couple 









msec 
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bubble , 
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71 

1380 

5 

1 

.5 


L 

ol 

6.6 
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44.5 

1 

59 
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11 

1 

.5 
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.1 

6.6 

16.8 

13.8 

2 

64 

2160 

1.4 

1 
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L 

.1 


6.6 

16.8 

15.5 

3 

17 

1020 

3 

1 

.1 


L 

.1 


6.6 

9.8 

8.8 

1 

23 

2100 

8 

1 

.5 


L 

.1 

0 

6.9 

2^9 

1 

36 

1740 

8 

1 

.1 


L 

.1 

0 

17.1 

15.3 

2 

99 

2190 

6 


.8 

1 

.1 


3 

46.6 

50.0 

3 
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1873 

3 

0 

.8 
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1 
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3 
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1 
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10 

1 

.5 
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1 
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10 

1 
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1 
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2 
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2 
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2 
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2 
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2 

0 

.8 
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3 
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1 

0 
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13 


.8 
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3 
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2 

60 
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3 
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.8 

0 


0 
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1 
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1 








109.3 
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2 
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2 
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117 
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13.2 
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a Bubble occurred after reference timer mark rather than before timer mark. 
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A voltage divider placed the voltage impressed on the timing lamp across 
a galvanometer. The lag of the square-wave timing-lamp ciruit in recording on 
the galvanometer was considered small. A constant five -frame timer-lamp 
advance with respect to the bubble history is caused by the position of the 
lamp in the camera (fig. 2). 

On the basis of these time relations, the reference time between optical - 
and oscillograph -recording systems may be summarized as shown in figure 5. If 
the reference time based on the optical system (t ro ,2 or t rQ i) equals the 

time measured on the oscillograph paper (t m 2 or ^rm l)j the a °tion occurred 
exactly as predicted from the observation 

tro,2 = X + ^ ~- (l) * * * 5 ’ + to + tm + cp (fig. 5(b)) 

tro,l =| ~ X + C g | ~ ^ “ tcp " <P (fig- 5(c)) 

where x is the number of frames from the reference timer mark to the initia- 
tion of the bubble, y is the number of frames it takes the bubble front to 
cross over the thermocouple location, cs is the camera speed in frames per 
second j and to and t<p are system response time and phase lag, respectively) 
cp is a 135° phase lag, which occurred in one run in the 60-cycle timer because 
of a nongrounded transformer. A segment of an actual oscillograph trace 
(fig. 6) is marked to illustrate these time delays. The measurements taken 
appear in table I. 

DISCUSSION 

When high-speed motion pictures of nucleate boiling (heat flux, approx. 
25,000 Btu/ (hr) ( sq ft)) were matched to the surface temperature profiles, a 

rapid temperature decrease assoc- 
iated with bubble growth was 
observed. This temperature de- 
crease, comparable to that 
noticed by Moore and Mesler 
(ref. 10), had four important 
characteristics : 

(l) The time relations in 

table I show that the drop occur- 

red as soon as the perimeter of 

the bubble base passed over the 
thermocouple junction. The 

error in these measurements is 
approximately ±2 milliseconds, 
a short time compared with the 
bubble life. 
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(2) The temperature drop indicates that the surface heat flux during 
growth is approximately 20 times the average heat flux electrically generated 
in the surface. These data were taken with the 0.0056 -inch -thick nickel-chrome 
alloy ribbon. 

(3) The temperature drop lasts for about 12 milliseconds, a time that is 
long in comparison to the heater relaxation period of 0.5 millisecond; 
therefore, the surface cooling mechanism beneath the bubble is considered to be 
continuous, with a duration of several milliseconds for the bubbles studied. 
These data were taken with the 0.001-inch-thick nickel-chrome alloy ribbon. 

( 4 ) The temperature reaches a minimum as the bubble approaches its 
maximum diameter. In all cases the temperature drop is completed before 
bubble collapse. 

The first and second observations of the temperature drop are consistent 
with an evaporative phenomenon only. The fact that such rapid cooling occurs 
directly beneath the bubble cannot be explained by forced convection. The 
third observation indicates that evaporation occurs continuously from a micro- 
layer as described by Snyder (ref. 5) and Edwards (ref. 6), rather than solely 
from the perimeter of the bubble base. Since evaporation in this experiment 
was limited by the thermal capacity of the heater ribbon, it is possible that 
under different surface conditions a dry spot might occur within the bubble 
base. The fourth observation demonstrates that surface cooling in no way 
depends on a quenching process; this conclusion, however, may not be valid 
below the subcooling range studied (5° to 20° F). Detailed temperature and 
heat-flux profiles of two typical bubbles are presented in figure 7. 

It can be concluded that the temperature drop is caused by liquid evapora- 
tion. In view of the experiments of reference 8, however, which showed that 
streams of air bubbles can raise the heat -transfer coefficient considerably, it 
is necessary to examine more closely the effects of bubble agitation. When a 
bubble leaves the heating surface, turbulence would be expected to be at a 
peak; yet, surprisingly, the heat flux is quite low and occasionally becomes 
negative. The reason for this adverse effect of agitation is that the surface, 
which has been cooled almost to T ga ^ by evaporation, is actually colder than 

the liquid that quenches it. Conversely, in the case of air bubbles, where 
there is no evaporation on the heater, quenching will take place over a hot 
surface; thus, turbulence will be a much more effective means of heat transfer 
than in actual boiling. This interpretation fits nicely with the results of 
reference 8, since the air bubbles did affect the heat flux strongly but less 
strongly than vigorous nucleate boiling did. 

Still, the evidence is not quite conclusive that the quenching mechanism 
is ineffective in all cases of boiling. The heater ribbon used in the present 
tests had a very small heat capacity, and during bubble formation all the 
available heat was withdrawn from the surface much more rapidly than it could 
be restored by electric generation. When boiling takes place on a more mas- 
sive solid, which has a higher thermal capacity and conductivity, it is pos- 
sible that the surface temperature after complete evaporation remains appre- 
ciably higher than that of the quenching liquid. Under these conditions and in 
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the case of a high degree of subcooling a quenching mechanism might become 
significant. Whether such rapid heat replenishment into the surface ever 
occurs in boiling remains to be proved. 

Even if the quenching mechanism is inoperable, microconvection in the 
liquid around the bubble could be an important path of heat transfer. As 
mentioned previously, it was demonstrated in the work of reference 3 that the 
thermal layer is greatly diminished over a distance of one bubble diameter 
around the nucleation site. In order to investigate the direct cooling effect 
of microconvection, the temperature profiles measured at three adjacent 
positions in the area of influence were compared (fig. 8 (a), (b), and (c)). 

In figure 8(a) the thermocouple is slightly less than one bubble diameter away 
from the nucleation site. No appreciable effect on heat transfer is seen dur- 
ing growth, and there is only a small rise in heat flux during collapse. In- 
creased pool turbulence probably accounts for this rise. In figure 8(b) the 
bubble base comes very near, but does not cover, the thermocouple junction; 
the results are similar to those of figure 8(a). In figure 8(c) the edge of 
the bubble base does pass over the thermocouple, and the rapid temperature 
drop associated with evaporation can be observed. Figure 8 demonstrates 
adequately that the direct cooling effect of microconvection due to bubble 
growth and collapse is small. Increased pool turbulence, by diminishing the 
thermal layer around the bubble , probably has a sizable indirect effect on the 
heat transfer. As pointed out in reference 3, the thickness of the thermal 
layer exercises extensive control over the nucleation capabilities of the 
surface . 

In terms of heat-transfer models, an evaporating microlayer at the base 
of the bubble is evidently responsible for most of the increased heat flux due 
to ebullition. 2 Neither microconvection nor surface quenching seem to be very 
important, although increased pool turbulence may well affect both surface 
cooling and nucleation appreciably. The contribution of ebullition to the to- 
tal heat transfer seems to depend mostly on the time-average fraction of the 
surface covered by an evaporating microlayer (i.e., bubbles in the growth 
stage). When the waiting period is many times the bubble life, free convec- 
tion will be the dominant heat-transfer mechanism. If, however, the waiting 
period equals the growth period, practically all the heat from the surface 
beneath a bubble column will be removed as latent heat. In this case, much of 
the heat supplied during the waiting period will not leave the solid but must 
be used to reheat the surface to a temperature sufficient for nucleation. 

As an illustration of where the evaporative mechanism becomes effective, 
consider the boiling case when an average of one-tenth of the surface area is 
periodically covered by bubbles. If the heat flux through the bubble base is 


3 An interesting calculation to approximate the thickness of the microlayer 
can be made by integrating the heat flux over the bubble residence time; enough 
heat is transfered during the temperature drop to evaporate completely a film 
8000 A thick at the bubble base . The actual layer is probably thicker than 
this because evaporation into the observed bubble was limited, not by the 
liquid available in the film, but by the small heat capacity of the heater 
ribbon. 
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generally 10 times the gross surface flux, evaporation and convection should 
have approximately equal effects on the heat transfer. When the surface 
becomes densely populated with bubble columns, it seems quite likely that 
latent heat transports as much as 70 to 90 percent of the total heat flux. It 
is reported in reference 12 that near burnout the fraction of the surface 
covered by bubbles was 0.4 or greater for distilled water; reference 10 also 
gives data in the 70 to 90 percent range that were calculated from heat flux. 


CONCLUSIONS 

The initiation of cooling due to bubble growth on a heating surface was 
investigated. It may be concluded from comparison of photographic and surface - 
temperature data that 

1. A very rapid temperature drop occurs when the perimeter of the bubble 
base passes over the thermocouple junction. 

2. A twentyfold increase over the average heat flux is observed directly 
beneath the bubble during the growth period; this surge in heat flux is 
completed as the bubble reaches its maximum diameter. 

3. There is an appreciable decrease in the local heat flux beneath the 
bubble during collapse. 

4. A small, but measurable, rise in heat flux occurs in the area of 
influence during bubble collapse. 

These results show that an evaporative heat -transport mechanism is much 
more important in nucleate boiling than has been indicated by previous experi- 
ments . In particular, strong support is given to the evaporating microlayer 
hypothesis of Snyder and Edwards. The microconvection and vapor-liquid dis- 
placement mechanisms of bubble agitation seem to have very little effect on 
heat transfer; this conclusion is independent of the heat flux though possibly 
not of extensive subcooling. Increased turbulence due to bubble collapse has 
a small direct effect on the heat flux, which may be appreciable in slow 
boiling. The contribution of bubble formation to the total heat transfer must 
be considered in relation to the entire ebullition cycle and will depend 
strongly on the time -average fraction of the surface covered by growing bub- 
bles . 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, January 29, 1964 


REFERENCES 


1. Rohsenow, W. M., and Clark, J. 
Transfer. Trans. ASME, vol. 


A. : A Study of the Mechanism of Boiling Heat 

73, 1951, pp. 609-620. 


17 


l 



2. Gunther, F. C., and Kreith, F.: Photographic Study of Bubble Formation in 

Heat Transfer to Sub cooled Water. Prog. Rep. 4-120, Jet Prop. Lab., 

1950. 

3. Hsu, Y. Y., and Graham, R. W.: An Analytical and Experimental Study of the 

Thermal Boundary Layer and Ebullition Cycle in Nuclear Boiling. NASA 
TPT D-594, 1961. 

4. Forster, K* E., and Grief, R.: Heat Transfer to a Boiling Liquid; Mecha- 

nism and Correlations. ASME Paper 58-HT-ll, 1958. 

5. Snyder, N. R., and Edwards, D. K.: Post Conference Comments; Summary of 

Conference on Bubble Dynamics and Boiling Heat Transfer. Memo. 20-137, 
Jet Prop. Lab., June 14-15, 1956, p. 38. 

6. Edwards, D. K.: The Role of Interphase Mass Transfer in the Mechanism of 

Nucleate Boiling. M. S. Thesis, Univ. Calif., 1956. 

7. Bahkoff, S. G., and Mason, J. P.: Heat Transfer from the Surface of a 

Steam Bubble in Turbulent Subcooled Liquid Stream. AIChE J., vol. 8, 
1962, pp. 30-33. 

8. Mixon, F. 0., Jr., Chon, W. C., and Beatty, K. 0., Jr.: The Effect of 

Electrolytic Gas Evolution on Heat Transfer. On: Heat Transfer. 

AIChE Chemical Engineering Progress Symposium Series, vol. 56, no. 30, 
1960, pp, 75-81. 

9. Westwater, J. W.: Boiling of Liquids. Advances in Chem. Eng., vol. I., 

Academic Press, Inc., 1956, pp. 1-76. 

10. Moore, F. D., and Mesler, R. B.: The Measurement of Rapid Surface Tem- 

perature Fluctuations During Nucleate Boiling of Water. AIChE J., 
vol. 7, 1961, pp. 620-624. 

11. Carslaw, Horatio S., and Jaeger, John C.: Conduction of Heat in Solids. 

2nd. ed., Univ. Press (Oxford), 1959, pp. 100-102. 

12. Donald, M. B., and Haslam, F.: Mechanism of the Transition from Nucleate 

to Film Boiling. Chem. Engr. Sci., vol. 8, 1958, pp. 287-294. 


18 


NASA -Langley, 1964 E-1432 






“The National Aeronautics and Space Administration . . . shall . . . 
provide for the ividest practical appropriate dissemination of information 
concerning its activities and the results thereof . . . objectives being the 
expansion of human knowledge of phenomena in the atmosphere and space. 

— National Aeronautics and Space Act of 195B 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles or meeting papers. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



